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Abstract 



All ALEPH measurements of branching ratios of r decays involving kaons are 
summarized including a combination of results obtained with Kg and K® detec- 
tion. The decay dynamics are studied, leading to the determination of contribu- 
tions from vector K*(892) and if* (1410), and axial-vector ifi(1270) and ifi(1400) 
resonances. Agreement with isospin symmetry is observed among the different final 
states. Under the hypothesis of the conserved vector current, the spectral function 
for the KKir mode is compared with the corresponding cross section for low en- 
ergy e + e~ annihilation, yielding an axial-vector fraction of (94^|)% for this mode. 
The branching ratio for r decay into all strange final states is determined to be 
B(t~ ->• X~(S = -l)v T ) = (28.7 ± 1.2) x 10~ 3 . The measured mass spectra of 
the strange r decay modes are exploited to derive the S = — 1 spectral function. 
A combination of strange and nonstrange spectral functions is used to determine 
the strange quark mass and nonperturbative contributions to the strange hadronic 
width. A method is developed to avoid the bad convergence of the spin zero hadronic 
component, with the result m s (M%) = (176 1 57 ) MeV/c 2 . The evolution down to 
1 GeV gives m s (l GeV 2 ) = (234 1^) MeV/c 2 . 
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1 Introduction 



Hadronic r decays provide a clean environment for studying the physics of hadrons, which 
are produced via W exchange, i.e., from the QCD vacuum. For this reason, the global 
properties of hadronic systems in r decay are described by fundamental quantities, called 
spectral functions, which measure the transition probability to create hadrons out of the 
vacuum as a function of the hadronic mass. Due to their unambiguous theoretical and 
experimental definition, r spectral functions provide information on hadron dynamics in 
an interesting mass region which is dominated by resonances and leads to the simpler 
asymptotic QCD regime. The ALEPH Collaboration has already published analyses 
of the nonstrange vector [ l| and axial- vector [ § spectral functions. The nonstrange 
vector spectral function can be compared to the corresponding cross sections in e + e~ 
annihilation [ p]] in order to test isospin invariance of the electroweak vector current, 
formerly called the conserved vector current (CVC) hypothesis. Detailed QCD studies 
have been performed using the nonstrange vector and axial-vector spectral functions, 
resulting in a precise determination of the strong coupling constant a s (M^) [ In this 
paper, attention is turned to the r decays involving kaons, which are of importance in 
order to address issues in both strange and nonstrange decay dynamics [ [| |J . 

Recently, ALEPH completed the measurement of r decay branching ratios involving 
kaons [ |5j, || [7], |J, as summarized in Table [l], which allows a comprehensive study of 
the strange sector and of some aspects of the dynamics in the final states with a KK 
pair. These analyses involve the detection of charged kaons and neutral kaons (K®, 
Kg). Assuming CP invariance the Kg and K\ branching ratios can be averaged. The 
corresponding results are given in the right-hand column of Table [l|. Although the 
strange final states are Cabibbo suppressed, the precision achieved in these measurements 
enables many significant tests, ranging from p — r lepton universality to hadron dynamics, 
resonance production and QCD analyses. In addition, the results are checked for 
consistency with isospin relations between the rates observed for different final states. 

One of the outstanding issues in the nonstrange sector is the vector and axial-vector 
composition of the KKir final state. This was a limitation in the QCD analysis of Ref. [ 
New information on this point is obtained from the ALEPH measurement of all 
relevant decay modesQ K^K^ir', K s K s ix~, K + K~-k~ and K°K~7r°. In addition, the 
availability of the corresponding isovector cross section measured in e + e~ annihilation 
and assuming CVC allows an independent determination of the vector fraction in the 
(K Ktt)~ v t channel, with a significant improvement in accuracy. 

Because G-parity cannot be defined for the strange hadronic states, it is difficult 
to separate vector and axial-vector contributions. Furthermore, the relatively low 
statistics do not permit different spin-parity states to be unfolded from the overall decay 
distributions. Experimentally, the Ktttt system is the most complex one to deal with, 
since it contains contributions from two axial- vector mesons, i^i(1270) and i^i(1400), 
and one vector meson, K*(1A10). However, the K*(1A10) state can be measured in the 
(K-k)~v t channel by its interference with the dominant K*(892) production, analogously 
to the p(770) - p(1450) situation in the tttc final states [ |. Unlike the K^UOO) and 

throughout this paper, charge conjugate states are implied. 
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Table 1: Summary of branching ratios for t decays involving kaons from ALEPH data. 
Channels with neutral kaons are measured separately using their K\ and K s components. 
Modes are split into the strange (S = —1) and nonstrange (S = 0) sectors of r decays. 
Due to the limitations of statistics and particle identification, the net strangeness of the 
K°h~h + h~ mode has not been determined (h stands for n or K). The last column presents 
the average of K s and K\ results for final states with K° and K~K°, and the sum of all 
components for final states with K°K°. 
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If* (1410), the axial- vector K i (1270) decays significantly into Kp. And so fits to the 
invariant tttt mass distributions [ [|, ^, [/], which measure the ifp fraction, determine 
the total if i (1270) contribution. A fit to the invariant Ktttt mass spectrum provides 
additional information on the resonance content. It is important to separate the vector 
and axial-vector contributions in the strange sector of r decays to obtain information on 
the size of the nonperturbative QCD part in the r hadronic width. This has been observed 
to be very small in the nonstrange case [ |2| . 

Independently of the resonance structure, the total strange spectral function is 
determined, without separating vector and axial-vector contributions. Similar to the 
nonstrange sector, the r strange spectral functions are key ingredients for QCD studies 
and provide the possibility to follow the results as a function of a variable "r mass" [ . 

One of the free parameters of the Standard Model, the strange quark mass m s , appears 
in many phenomenological calculations, such as in the prediction of the CP-violating 
kaon parameters e'/e [ [|. The strange r decay rate is sensitive to m s [ [l(| which can be 
experimentally determined. A breaking of chiral symmetry, induced by the relatively large 
strange quark mass, introduces a mass dependence into the perturbative QCD prediction 
of the total strange hadronic width of the r. Using the total rates and moments of the 
respective spectral functions, a combination of strange and nonstrange modes can be found 
which cancels the dominant massless perturbative contribution, enabling the strange 
quark mass to be derived by means of a combined fit with nonperturbative contributions. 
An error in the theoretical calculation [ |TTJ of the second-order perturbative term was 



recently pointed out [12]. After correction, a problematic convergence behaviour of the 



perturbative series is observed [ [12], |13|, |T4}| , leading to a larger theoretical uncertainty 



on m s than had previously been assumed and to a bias in the central value obtained [ 
[3|] . In this analysis, particular attention is devoted to this question and a new method is 
introduced [ [TH] to avoid the convergence problem. 



2 \jl — t lepton universality 

Lepton universality requires that the W~ boson coupling to any lepton pair lv\ be of the 
same strength gf, this hypothesis is in agreement with studies in the leptonic and hadronic 



sectors [ |16[ . Lepton universality can be tested in the strange sector by comparing the 
decay widths of r~ — > K~v T and K~ — ► pTv^. The following ratio, in which the common 
CKM matrix element and the decay constant fx both cancel, is precisely predicted by 
theory: 



Rt/K = 



T(t~ -> K~v r ) 

T(K~ -> /i"^) 

B(t- -> K-u t )t k 

B{K- -> p~u^)T T 

Igl Ml (1-M 2 K /M 2 T )\ 



1 + 6R T/K ) 



2glM K Ml{l-Ml/Miy 

with the radiative correction 6R t /k = (0.90 ± 0.22)% [ 0, B(K~ -> fTu T ) = B K = 
(63.51 ± 0.18)%, r K = (1.2386 ± 0.0024) x 10~ 8 s, and r T = (290 ± 1.2) x 10~ 15 s [ p 
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Assuming ji — r universality (g T = g^), Eq. ([I]) predicts 

B(t~ -> K-i/ r ) = (7.14 ± 0.02 t/l ± 0.02 Ba . ± 0.03 r J x 10~ 3 , (2) 

where only sizeable contributions to the errors are given (the theoretical uncertainty is 
from radiative corrections). 

The comparison of (0) with the ALEPH measurement 

B(t~ -> K~v T ) = (6.96 ± 0.25 ± 0.14) x 10~ 3 (3) 

yields 

g r /9v = 0.987 ±0.021, (4) 

which agrees with the universality assumption within one standard deviation. 

The measurement (H) can also be exploited to determine fx if lepton universality is 
assumed to hold. One obtains^ 



where G F = (1.16639 ± 0.00002) x 10~ 5 GeV" 2 [ [18| is the Fermi coupling constant, and 
radiative corrections are ignored. This results in 

f K = (111.5 ±2.3 ±0.9) MeV , (6) 

where the first error is from the uncertainty on B(t~ — ► K~v T ) and the second from I4 S . 



3 Study of resonance structure and dynamics 

The decay dynamics of r decays involving kaons are investigated via the invariant mass 
spectra. Due to the small statistics, decay modes other than (Kir)~, (if7r7r)~ and (KKtt)~ 
are not discussed here. Even for these, only a qualitative approach based on simple 
descriptions can be considered. 



3.1 Mass spectrum for r — ► (Kir) v r 

Three measured modes, K~tt°, Kgir~ and K\-k~ , are used to study the invariant mass of 
the K7i system. The if* (892) ~ dominance is well established in all the cases [0, 0]. For 
the first two modes, where the mass resolution is good enough, some excess is observed at 
large mass over and above the expected if* (892) ~ tail. This could arise from a possible 
if*(1410) vector state contribution, expected to be dominated by if* (892) — if* (1410) 
interference as already observed for p(770)— p(1450) in the tc~7t° system [ |I|. Theoretically, 
the decay rate for r~ — > (Ktt)~u t relative to the electronic width is expressed as [ [19], ^(J 




x(s ~ Mlf (s - Miy 2 s-V(s)f , (7) 

2 The decay constant fp should be scaled by a factor of \/2 to compare with the value given in Ref. [ 
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where M± = M K ± M v and the form factor includes both if* (892) and if* (1410) 
contributions, namely 



f(s) = T —[BW K .( f m)(s) + PBW K * (im (s)} , 



where (3 is assumed to be real [ ^(J . The two Breit-Wigner propagators have the form 



/If 2 



with an energy dependent width T(s) 




p(s) = ^(s-Ml)(s-Ml) , (11) 

where I is the K — n angular momentum, i.e., I = 1 for p wave. 

The if* (1410) search is performed using only the K~n channel as this has the lowest 
background in the interference region and the best mass resolution. Figure [l] shows the 
mass spectrum after background subtraction and bin by bin correction for both resolution 
and acceptance. The fits using the above parametrization are listed in Table [|. If the 
mass and width of the if* (892) are left free (Fit 1), the fitted values are found to be 



consistent with the world average [ i~8fl . With the fixed K* parameters, the final value is 



obtained (Fit 2) using the K 7r° data, where the systematic uncertainty from the non-iT7r 



events is negligible. A fit to the combined K 7r° and K^tc invariant mass distributions 



is also given (Fit 3). 

The measurement of /3 — —0.11 ± 0.05 is, as expected from SU(3) flavour symmetry, 
in agreement with the p(1450) contribution ((3 = —0.094 ± 0.007 from the ALEPH 
measurement of the decay r~ — > 7T~7T v T [ [I]]). Low statistics do not allow much more 
information to be derived from the fit; however, leaving the phase of (3 as a free parameter 
yields a value (180 ± 55)°, consistent with the assumption made. 

The relative K"*(892) and K*(1410) "diagonal" contributions to the t~ -> (Ktt)-u t 
branching ratio determined from the value of (3 are found to be (97.5l{;4)% and (0.8^q;5)%, 
respectively, the interference contributing (1.7^1)% for the remainder. According to 
Ref. [ [18] the decay fraction of fT*(1410) into Kn is (6.6 ± 1.3)%, and for this analysis the 
rest is assumed to decay into Ktxtx final states. After correction for the r kinematic factor 
(see Eq. (|34"D), this fraction becomes (7.2 ± 1.4)% for K*(1410) produced in r decays and 
leads to 

B(t~ -> K*(UW)~v T ) = (1.5t\i) x 10 " 3 > ( 12 ) 

using the Kir branching fraction from Table |I|. 

In this analysis, a contribution from the scalar Kq (1430) has been ignored. No evidence 
for this state is seen in the Kir mass plot. Nevertheless, a fit with this contribution was 



performed using the resonance parameters from Ref. [ . A 95% confidence level upper 
limit of 0.5 x 10~ 3 is derived for the corresponding branching ratio. 
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Figure 1: The invariant mass spectrum for t~ — > K~Tt°i> T after background subtraction, 
bin migration and acceptance corrections. The fit is described by the solid curve, taking 
into account both if* (892) and if* (1410) contributions with interference. The dashed and 
dotted curves show the if* (892) and if* (1410) parts. 



Parameter 


Fit 1 {K-7i°) 


Fit 2 (K-7T ) 


Fit 3 (K-tt° + K° s 7i-) 


M K * {892) (GeV/c 2 ) 
rV (892 ) (GeV/c 2 ) 

P 

M K * {im (GeV/c 2 ) 
r^ ( i4io) (GeV/c 2 ) 


0.895 ± 0.002 
0.055 ±0.008 
-0.06 ±0.07 

1.412 

0.227 


0.892 
0.050 
-0.11 ±0.05 
1.412 
0.227 


0.892 
0.050 
-0.12 ±0.04 
1.412 
0.227 


X 2 /ndf 


13.0/19 


16.0/21 


19.8/21 



Table 2: Fit results for r~ — > K^-k q v t and t~ — > Kgir~v T decays. The errors given are 
statistical only. Because the background of r~ — > K~K°i> T dominates over the higher 
K s ir~ invariant mass tail, the background subtraction relies on the model prediction [\2f\j 
and introduces additional systematic uncertainties in the K s n~ system. The fixed values 
for the masses and widths of the if* resonances are taken from Ref. [ TUrf. 
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3.2 Mass spectra for r — > (Kitty) v t 

The (fT7r7r)~ decay occurs in the four final states i^~7r + 7r~, K^tt^tt , K®tt~tt° and 
fC _ 7r°7r . Its invariant mass spectrum as well as the spectra for the Ktt and tttt subsystems 
reveal details of the underlying decay dynamics. 



3.2.1 Contribution from K*(U10) 

The fitted contribution of if* (1410) found in the previous section has consequences for 
the (Knrty final state. The fraction (92.8 ± 1.4)% of if* (1410) decaying into Ktxtx is 
used to estimate the vector Kttti component, while the remainder is attributed to the 
axial- vector contribution dominated by i^i(1270) and i^i(1400). For the vector current 
part, one obtains 

B{ T - -> K*{im)-u T -> (Kttk)-v t ) = (1.4±5;1 ±g:2) x 10" 3 , (13) 

where the first uncertainty comes from the fit to the Ktt invariant mass, while the 
second uncertainty arises from the possibility for the X*(1410) to decay into Krj. After 
subtracting this vector component, the axial- vector contribution to the Ktttt decay modes 
is 

B A {r- -> (Ktt7t)-p t ) = (4.61^) x 10~ 3 . (14) 



3.2.2 Kp and K tt fractions 

The Ktt and tttt mass spectra are used to search for the intermediate states K tt and Kp. 
The Kp fraction can reveal the presence of the axial vector ^(1270) resonance which 
decays about 50% of the time to Kp [ [L8|; unlike the ifi(1400) which decays almost purely 



to K tt. For practical reasons related to the shape of the different backgrounds, the most 
efficient way for separating the K tt and Kp intermediate states is to fit the tttt invariant 
mass spectra where background and non- Kp contributions are essentially located at lower 
masses. 

No excess in the p mass region is found in the K~tt°tt° mode, as expected, while a p 
signal is observed in the K~tt + tt~ mode [ |J and also in both the K^tt^tt [ |6| and K1tt~tt° 
modes [ |7J . With the assumption that the Ktttt final states proceed only via an incoherent 
superposition of the intermediate states Kp and K it, a sum of a p Breit-Wigner signal 
with the r kinematic factor and the shape of the K tt reflections obtained from simulation 
is used to fit the tttt mass spectrum after subtracting the tioyi-Ktttt background. It yields 
the fractions f K - p0 = (39 ± 14)% for the K-n+n- mode0 and f w = (74 ± 13)% for 

the K°tt~tt° mode. The sum of the fits to the tttt invariant mass spectra is shown in 
Fig. 0. These results indicate that the {Kp)~ intermediate state plays an important 
role in the Ktttt mode. The smaller value in the K~tt + tt~ mode is understood from 
isospin considerations as discussed in Section 5. The sum of decays involved in the Kp 
intermediate state gives B(t~ — > (Kp)~u T ) = (3.25 ± 0.67) x 10" 3 . 

3 This value differs slightly from the published result (35 ± 11) in Ref. [ |5| which was obtained with a 
p lineshape inadvertently not corrected for the r kinematic factor. 
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Figure 2: The tttt invariant mass spectrum for t~ — > (Ktt7t)~v t after background 
subtraction. The overall fit, the expected p signal and K*tt reflection are shown in solid, 
dotted and dashed histograms. 



The branching fraction of the if i (1270) into Kp is taken from Ref. [ y| and becomes 
(52 ± 7)%, when taking into account the r kinematic factor. A first estimate of the 
branching ratio to the if i (1270) is therefore derived to be 

B{t~ -> tff(1270)i/ T ) = (6.3 ± 1.5) x KT 3 (Kp fractions) , (15) 

where all the Kp states are assumed to originate from Z£ 1 (1270) decay only. 



3.2.3 Fit to Kittt resonances 

The description of the Kittt mass spectrum in terms of the resonant states ifi(1270), 
ifi(1400) and if* (1410) is rather complicated because several intermediate states are 
involved, as discussed in Ref. [ ^2). Therefore, the shape of the respective mass 
distributions is obtained from the Monte Carlo simulation, taking into account all the 
known intermediate states. The resonance parameters for the K i (1270), if i (1400) and 
if*(1410) are fixed at 1.273 GeV/c 2 , 1.402 GeV/c 2 and 1.412 GeV/c 2 for the masses and 
at 0.15 GeV/c 2 , 0.174 GeVc 2 and 0.227 GeV/c 2 for the widths. The chosen width for 



the if i (1270) is different from the value in Ref. [ IS| , because of the parametrization 
adopted: it should be interpreted as an effective width, resulting from the opening of the 
different decay channels. It has been checked that the resulting line shape agrees well 



with the existing data [ 22fl . However as the value for the width can have a strong 
effect on the relative fractions in the If i (1270) decay, an uncertainty of ±0.05 GeV/c 2 is 
assigned to it. The mass resolution and the statistics are not sufficient to separate the 
if i (1400) and if* (1410) states, and so an effective resonance is used instead, averaging 
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the parameters of the two states. To obtain the shapes including all the intermediate 
states, the 1^(1270) decays into Kp } ir*(892)vr, ir*(1430)vr and iT/ (1370) are generated 



separately and then combined according to the relative fractions given in Ref . [ . Since 
the r kinematic factor modifies the branching ratios, the fractions for the different final 
states are recomputed for ^(1270) production in r decays, yielding (52 ±7)%, (25 ±8)%, 
(8 ± 1)%, (14 ± 3)% and (1 ± 1)% for the Kp, K*(892)tt, K*(U30)n Kuo and Kf final 
states, respectively The fit is shown in Fig. [|, giving y 2 /ndf = 14.1/22 and a fraction 

/x l( i270) = (41 ±19±15)%, (16) 

where the first error is of statistical origin and the second comes from the i^i(1270) 
width. Because the total branching ratio for the Ktctt modes is B(t~ — > (K7T7t)~v t ) = 
(5.97 ± 0.72) x 1(T 3 , the expected branching ratio for i^i(1270) is 

B(t~ -> K{(1270)u T ) = (2.9 ± 1.7) x 10~ 3 (Kirir fit) , (17) 

taking into account the fact that (86 ±3)% of the i^i(1270)'s decay into Kirn final states. 



The values (flq) and (|T7| ) obtained from two independent estimates are in fair 



agreement and their combination yields 

B(t~ -> Ki (1270)i/ T ) = (4.8 ± 1.1) x 10~ 3 (combined) . (18) 

From this value, the r decay into K~uj through iff (1270) is estimated (using the 
branching ratio in Ref [ |IB|| ) to be (0.67 ± 0.21) x 10~ 3 , indicating an axial-vector 
contribution of Ba{t~ — ► K~7T + 7T~7r°v T ) = (0.59 ± 0.19) x 10~ 3 in the measured Ktttttt 
decay modes. Comparing the results in Eq. ([14]) and Eq. (|l~8|) , one obtains a branching 
ratio B(t~ — > (1400)z/ r ) = (0.5±1.7) x 10^ 3 , which is somewhat reduced in comparison 
with that of i^i(1270). 

3.3 Mass spectra for r — > (KKn) ~v T 

Almost all possible final states for (KKn)~ are studied by ALEPH (Table |1]): K + K~tt~ , 
K° s K° s tt-, K° s K° l tt-, K~K° s tt° and K-K° l tt°. The K+n' mass plot in the K + K~n~ mode 
shows a clear K*(892)° signal, and a fit yields a fraction of (87 ± 13)% for the K*K~ 
component [ A strong i^*(892)~ signal is also found in the K~K° L ir° mode with an 
estimated lower limit of 86% at 90% confidence level for the K (K*) component [ [/J, 
confirming the conclusion of K* dominance in the KKir mode. 

Both the K + K~tc- and K~K^tt° modes are used to measure the total KKtc invariant 
mass distribution, scaling their rates by a factor of two. Figure ^ shows the combined 
KKtt invariant mass-squared distribution after correcting for resolution and efficiency. 



The Monte Carlo expectation [ |2"3"|, based on ai dominance, is observed to be consistent 
with the data. 



4 Tests of isospin invariance 

Isospin relations relevant to r decays into kaons have been established in Ref. [ ^4], |25], |2f 
The set of results obtained by ALEPH allows for a complete test of these relations. Apart 
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M 2 (Ktut) (GeV/c 2 ) 2 



Figure 3: The invariant mass-squared distribution for r~ — > {K-k-k)^v t after efficiency 
corrections for the different final states and background subtraction. Data are shown by the 
dots with error bars. The fit is described by the solid curve, while the expected i^i(1270) 
and Zfi(1400) contributions are separately depicted by the dashed and the dotted curves. 
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Figure 4: T/ie invariant mass-squared distribution for t~ — > (KKn)~v T after efficiency 
corrections for the different final states and background subtraction. Data are shown by 
the dots with error bars. The histogram represents the model expectation f^SjJ. 
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Figure 5: The published branching ratios for r — > (if7r) v T WV , using all three 

measured modes and assuming I = 1/2. 

from testing isospin invariance in the hadronic states accessible in r decays, this procedure 
can be helpful as a consistency check among the different final states, regardless of the 
decay dynamics. 

4.1 The decay r~ — > (Kir)~v T 

The branching ratios of the three final states of the decay r~ — > (Kir)~v T , respectively, 
K~7T°, Kgir~ and K^-k~ , are presented in Table |1|. A study of the corresponding invariant 
mass spectra shows that the decay is dominated by if* (892), with a small component 
of if* (1410) visible through its interference with if* (892) as presented in Section 3.1. 
Isospin symmetry for an i = 1/2 state, together with the assumption of equal ifg and if£ 
contributions, constraints the rates for the three final states to be equal up to small mass 
corrections. Taking this into account, the above three measurements are well consistent 
with isospin symmetry, giving a x 2 /ndf = 0.9/2. Using isospin symmetry as a constraint 
yields the total branching ratio 

B(t~ -> (Ktt)~p t ) = (13.60 ± 0.62) x 10~ 3 . (19) 

A comparison between this analysis and the other measurements is shown in Fig. [5[ 

From the fit performed in Section 3.1, the branching ratio for if* (892) ~ production is 
found to be 

B(t~ -> if *(892)~z/ r ) = (13.26 ± 0.63) x 10" 3 . (20) 
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Figure 6: The decay fractions for K~tt + tx~v t and K°tt°it~i> t given by the standard error 
ellipse (39% probability). The isospin constraint is shown by the straight line. The 
triangles indicate possible intermediate states dominated by K tt, Kp or K~(tttt)i = o. 

4.2 The decay r~ — ► {Kitit)~v t 

For a total isospin I = 1/2 and the two possible isospin values 1^ = 0, 1 for the tttt system, 
a relation can be established among the branching ratios of the K°tt~tt°, K~tt + tt~ and 
K-tt°tt° modes [ gjj: 

B(t~ -»■ K-n + 7i-v T ) = 2 B ( T ~ K°K-n°v T ) + 2B(t~ -»■ K-tt°tt°u t ) . (21) 

In order to examine to what extent the above isospin constraint is satisfied by experiment, 
it is convenient to use the relative fractions of the total Ktttt branching ratio for each 
individual mode. This gives /(^vr-vr ) = 0.55 ± 0.06, f(K->K + iT-) = 0.36 ± 0.06, and 
f(K~7T°7T°) = 0.09 ± 0.04 from the ALEPH measurements with a correlation coefficient 
of —0.8 between f(K°7r~7r°) and f{K^-K + -K^). The corresponding error ellipse is given in 
Fig. H showing that the isospin relation (|2l"D is satisfied within one standard deviation. 

The isospin relations can be tested in more detail if the different intermediate states 
are separated. Three possible intermediate states occurring in the (Kttti)^ final state are 
considered: Kp and K~ (irn) I=0 corresponding to I n7r = 1 and 0, respectively, and K tt, 
with the expected relative rates in the different final states as given in Table |3]. It turns 
out that a pure K~(tttt)j =q is excluded by the measurement, while a combination of Kp 
and K n can fit the data well. This is in agreement with the findings of the resonance 
structure analysis performed in Section 3.2.2. A further check of the isospin symmetry 
is achieved taking advantage of the Kp and K it separation provided by the p resonance 
fits performed in the two relevant channels. The results given in Table [| are in agreement 
with the expectations. 
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Mode 


K 7T + 7T 






X 2 /ndf 


Kp 


1/3 


2/3 





7.0/3 


K*n 


4/9 


4/9 


1/9 


2.8/3 




2/3 





1/3 


88.7/3 



Table 3: The relative fraction for each Ktttt mode from isospin symmetry for the possible 
intermediate states. The test of these fractions is given by the y^jndf values for the 
hypothesis that the given intermediate state dominates. 



Mode 


K-n+Ti- (10~ 3 ) 


K°n-n° (10- 3 ) 


K'7r°n° (10- 3 ) 


X 2 /ndf 


Kp 

K*7l 


0.83 ±0.35 
1.31 ±0.41 


2.42 ±0.57 
0.85 ±0.45 




0.56 ±0.25 


0.7/3 
1.8/3 



Table 4: Test of isospin symmetry for different intermediate states. The separation of the 
branching ratios for Kp and K it is based on the measurements of the Kp fractions. The 
X 2 /ndf values represent the consistency with the expected relative ratios given in Table^. 



4.3 The decay r — ► (Ktttttt) v t 

The final states considered include K^-n^-n^-n^ , A'-7r°7r 7r°, K°h + h~h~ and K°tt~tt it , in 
which the K~r\ contributions are excluded. Since the net strangeness for the K°h + h~h~ 
modes is not determined, K°7r + 7T~7T~ dominance is assumed. Isospin invariance for the 



3n system leads to only three symmetry classes [ |2g, [28[: two with i 37r = 1, as shown 
in Table |j, and one with I 3n = 0, expected to be dominated by the K~oj state. An 
axial- vector contribution from the decay iff (1270) — > K~& — ► K~7t + 7t~7t , amounting to 
(0.59 ± 0.19) x 10~ 3 , was deduced in Section 4.2 and is subtracted from the K~7t + 7t~tt 
mode. 

The expected relative contributions from the 1^ = 1 classes are compared to data in 
Table |5|. The present experimental precision does not allow definite conclusions. However, 
the different isospin configurations can be used to derive a better estimate of the total 
rate. The two isospin possibilities lead to consistent values (see Table [5]) and the average 



using the Fermi statistical weights [ 28| yields (excluding the K u mode) 

B(t~ -> (Kirinr)-v T ) = (0.76 ± 0.44) x 10~ 3 . (22) 



4.4 The decay r — > (KKtt) v t 

For this mode, the decay dynamics are investigated using the branching ratios for all the 
relevant final states: K + K~tt~, K^K~tt°, K^K'tt , K° s K° l ti- and K° s K° s ir-, as shown 
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Class 


K 7T + 7T 7T° 


if-7r 7r°7r° 


K°7T + ir IT 


K°TT-TT 7r° 


5(Z3vr) (10~ 3 ) 


X 2 /ndf 


(300) 


2/15 


1/5 


8/15 


2/15 


0.65 ±0.39 


1.8/3 


(210) 


1/3 





1/3 


1/3 


0.83 ± 0.46 


3.0/3 



Table 5: The relative fractions for each Ktittti mode expected by isospin invariance for 
I 3n = 1. The three numbers {nin 2 n^) label Pais' isospin classes The test of isospin 

symmetry is given by the x 2 /ndf values for the hypothesis that the given intermediate 
state dominates. 



Mode 


K-K+7T- 


K° s K° L n- 


K° s K° s 7r- + K° L K° L 7T- 


K-K°7i° 


X 2 /ndf 


K*K 


1/3 


1/6 


1/6 


1/3 


2.1/3 




1/4 


1/4 





1/2 


16.0/3 


n(KK) I=0 


1/2 





1/2 





99.0/3 



Table 6: The relative fraction for each KKtc mode in view of isospin relations for three 
possible intermediate states. The test of these relations is given by the x 2 / n df values for 
the hypothesis that the given intermediate state dominates. 



in Table [!]. If the contribution of second-class currents is ignored, the rates for K + K it 
and K K°tt~ are expected to be equal [ |25|] , in agreement with the measurement 

*< T '-*^*>= 0.94 ±0.27. (23) 

The possible intermediate states appearing in these decay modes are considered to be 
K*K + K*K (denoted K*K), pit and 7r(KK)i = o, corresponding to the fractions shown in 
Table [!]. The test of isospin symmetry is consistent with the K*K intermediate state as 
illustrated in Fig. [7|, in agreement with the conclusion from the analysis of the Kir mass 
spectra. 

4.5 The decay r~ — > (K Ktttt)^u t 

The limited experimental information (see Table [D on KKnn final states makes it 
impossible to study this channel in any detail. Nevertheless, some conclusions can 
still be drawn. First, the rates for r~ — > K + K~7t~tt°v t and r~ — > K q K q -k^-k g v t are 
consistent within their uncertainties, as expected by isospin symmetry and the absence 



of second-class currents [ 25]. Then, the relative fractions observed in the different modes 



can be compared with the expectations from the possible isospin configurations [ [26 
with a strong preference for the = 0, I W7r = 1) states and to a lesser extent 

(Irk = 1> ^irir — 1)- The weighted average over the two non-zero modes yields an 
estimate for the total rate 

B(t~ -> {KKiwYvr) = (0.5 ± 0.2) x 10" 3 . (24) 
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Figure 7: T/ie decay fractions for K Q K~n Q v T and K~ K + Tt~i> T , given by a standard ellipse 
(39% probability), are compared to the isospin constraint on the KKir mode, shown by 
the straight line. The triangles indicate the intermediate states dominated by tt(KK)j = q, 
K*K or pix. 



5 Vector and axial-vector separation 

The separation of the r final states into vector and axial-vector currents provides 
important information for QCD studies. Any deviation between the inclusive sums of 
vector and axial- vector hadronic branching ratios is necessarily a result of nonperturbative 
phenomena. In this section, all the information obtained on the strange decay fractions 
concerning their V/A character are summarized in order to determine the corresponding 
inclusive rates. In addition, the separation of vector and axial-vector components is 
presented for the nonstrange KKn mode. 



5.1 r — > (Knit) v T 

The branching ratios for the strange sector of r decays are globally reconsidered in Table [7], 
with the respective current contributions unambiguously separated for the K~, (Kn)~, 
K~t] and (Knii)~ modes. In the (K3n)~ mode, the contribution expected from the 
axial- vector i£i(1270) decay into K~u is subtracted. The total (K4it)~ and (i^57r)~ 
contributions are estimated on the basis of the branching ratios for the 5n and Git modes [ 



18[ with Cabibbo and kinematic suppression, even though the resulting total contribution 
is small in comparison to the precision achieved in the measurement of the strange sector 
of r decays. For the non-K~uj part of (K3n)~ , the (KAtt)~ and the (K5n)~ contributions, 
it is assumed that vector and axial-vector contribute equally with 100% anti-correlated 
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MnHp 


h t , no- 3 ! 






1 1 




f. Qf. _i_ n 90, 


f. Qf. _i_ n on 




13 60 ± 62 




13 60 ± 62 


(K2n)- 


1.391^? 


4.58^:1 


5.97 ±0.73 


Kf(1270) -> if-w 




0.67 ±0.21 


0.67 ±0.21 




0.29^;1| 




0.29l°it 


(Z3tt)- 


0.38 ±0.53 


0.38 ±0.53 


0.76 ± 0.44 


(Z4tt)- 


0.17 ±0.37 


0.17 ±0.37 


0.34 ±0.34 


(K5ti)- 


0.03 ±0.10 


0.03 ±0.10 


0.06 ±0.06 


Sum 


15.86™ 


12.79 ItS 


28.65 ± 1.17 



Table 7: Branching ratios for vector and axial vector current contributions to the strange 
sector of r decays. The branching fractions for the KAm and K5tt modes are obtained 
from the measured branching ratios for the 57r and 67r in Ref. f\Wj] introducing the Cabibbo 
suppression and kinematic factors. Vector and axial vector currents in the K3n, KAtx and 
K5n are assumed to contribute equally. 



errors. The corresponding branching ratios for vector and axial-vector final states are 

B(t- -> V~(S = -l)u T ) = (15.9 tH) x 10~ 3 (25) 

and 

B(t- -> A-(S = -l)u T ) = (12.8 t{ A 7 ) x 10~ 3 . (26) 
The total branching ratio for all strange r decays 

B(t~ (V + A)-(S = -l)v T ) = (28.7 ± 1.2) x 10" 3 (27) 

is known with better precision. 

An estimate of the size of nonperturbative QCD contributions can be obtained from 
the ratio 

b(t- ^v~(s = -1)v t ) - b(t- ^a-(s = -ik) H 1 +10 



B(t- -> (V + A)-(S = -l)v T 



111?)%, (28) 



which is consistent with zero, but cannot rule out a significant nonperturbative 
component. In comparison, the value (1.7 ± 1.0)% was obtained in the nonstrange sector 
of r decays [ |2f . 

5.2 t~ (KKir)-v T 

In the analysis of the nonstrange spectral functions [ p], [J, fractions of (50 ± 50)% vector 
and axial-vector currents were assumed for lack of a better knowledge in r~ —>■ (KK-k)~p t 
decays. Due to the anti-correlation between vector and axial-vector, this uncertainty was 
important for the difference V — A. 
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In this analysis, two methods are considered in order to separate the vector and axial- 



vector contributions in the KKir mode. The first method based on the G-parity [ 
applies best to the cases of either a pure vector or a pure axial-vector state. The second 
method is to compare the inclusive V + A KKn spectral function of r decays with the 
corresponding isovector cross section measurements in e + e~ annihilation, providing the 
V part by means of CVC. 

The G-parity of the KKir system is determined by the KK system with G K -^ = 
( — l) Z A'7f . Since the K°K~tt° mode can only have I K -^ = 1, the corresponding 

decay width amounts to half of the total decay width for I K j^ = 1. The isospin 
configuration I K -^ = contributes to both the K°K°tt~ and K + K~n~ modes. In 
addition, the l K -^ values are always odd for K Q S K\ and even for KgKg or K\K\. The 



two following ratios are therefore defined [ |25 



and 



_ Y K%Kh - + T K o Kln - _ 2T° A + I* 



^K°K-TT° + Ty 



(30) 



where T VA and T VA are the decay widths for the l K -^ = and l K -^ = 1 states, respectively. 
For pure vector, one has 1/r = r , whereas r = tq for pure axial vector. The branching 
ratio measurements from ALEPH give 

r = 0.52l°;! (31) 

and 

r = 1.1918S - (32) 

Consequently, the probabilities for the KKn mode being a pure axial vector and a pure 
vector are computed to be about 25% and 50%, respectively. It is therefore not possible 
to draw a firm conclusion from this study. This is not too surprising as the KKn mode is 
observed to be dominated by K*K, in which case the previous test becomes degenerate. 

Another possibility is the comparison between the measured spectral function in r 
decay and the measured isovector cross sections in e + e~ — > KKtc. The application of 
CVC makes it possible to link the r decays into KKtt final states with e + e~ annihilation 
into isovector states KgK ± n T and K ± K T n°, via the relation 

i=i 4to 2 _ 

a e+e-^KK-K ~ g V l,KKnv T ' V 66 ) 

where a 1 ! 1 . ^ is the isovector cross section, a is the electromagnetic fine structure 

e+e~— >KKtt ° 

constant, and v 1 is the corresponding vector r spectral function. There is only 

isovector in r decay, while the isospin for the KKn system in e + e~ annihilation can be 
either or 1. By studying the Dalitz plot, the DM1 [ |29| and DM2 [ |30| experiments 



observed K*K dominance and concluded that the reaction e + e~ — > KKir is dominated 
by the isoscalar (resonant) amplitude, whereas a small isovector component is extracted 
from the interference of both amplitudes. 
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Figure 8: The V + A spectral function using the KKtt data from r decays compared to 
those derived from the isovector e + e~ cross sections. The e + e~ data are taken from DM1 [ 



H/ and DM2 [^Qj. 



The measured invariant mass-squared distribution of the KKtt system in r decay 
(Fig. [|) is transformed into the V+A spectral function (according to Eq. ( 531 ) in the next 
section). The cross section measurements for e + e~ — > KKtt [ |29], |3(| are also converted 
into a spectral function according to Eq. fl3"3"|). Since the DM1 and DM2 experiments 
only provide studies of K^K ± tt t , the full KKtt spectrum is obtained by scaling up the 
K$K tt^ contribution by a factor of three. Figure || shows the comparison between 
t decay and e + e~ annihilation. A fairly large excess is observed below 2.5 (GeV/c 2 ) 2 
compared to the 1 = 1 e + e~ data. The measured cross sections of e + e~ — > KKtt 
are translated by CVC into the expected vector branching ratio By(r~ — > (KKtt)~v t ), 
yielding (0.67 ± 0.27) x 10" 3 for the DM1 measurement and -(0.12 ± 0.26) x 10~ 3 for 
the DM2 measurement, which are averaged to (0.26 ± 0.39) x 10 -3 . This value is much 
smaller than the measured branching ratio B{t~ — > (K Ku)~ v r ) = (4.60 ± 0.50) x 10 -3 , 
yielding a dominant axial-vector fraction of (94i 8 )%. 

One can combine all the information sensitive to the V and A fractions: relative rates 
of different KKtt modes in r decay and comparison of total rates in r decays and 1 = 1 
e + e~ annihilation. It is found that the I K j^ = 1 component in the dominant axial-vector 



part is (75 ± 9)%, while in the much smaller vector part it amounts to (43 ± 28)%. These 
fractions agree with the expected value of 2/3 for K*K dominance, established for the 
axial- vector current in this paper, and for the vector current in e + e~ annihilation [ [29]. |3(J . 

A natural candidate to explain the dynamics of the KKtt mode is the decay 
di — > K*K. In fact the observed mass spectrum agrees well with the Monte Carlo 
prediction based on this model (Fig. |j), showing a sharp rise at the K*K threshold. 
Under the assumption of a\ dominance and using the branching ratio Bir — ► a\P T ) [ 
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one gets the branching fraction B(a\ — > K*K) = (2.6 ± 0.3)%. 
good agreement with the preliminary result from a partial wave analysis 
r 

width 



This value is in 
in the 
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im°ii v T channel, including the opening of the K*K decay channel in the a% total 



yielding B{a\ — > K*K) = (3.3 ± 0.5 ± 0.1)%. This observation provides an 
additional and independent argument for axial-vector dominance in the KKtt mode. 

The observed axial-vector dominance is in qualitative agreement with the models 
developed in Ref. [ pO, |23|] and in contradiction with those proposed in Ref. 
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6 The total strange spectral function 

Theoretically, the hadronic r decay width can be formulated in terms of the spectral 
functions [ |35| V\(s), ai(s) and ao(s) for the nonstrange part, and vf(s), af (s), Vq(s) 
and ciq{s) for the strange part, where s is the hadronic mass-squared. The notations 
v and a stand for vector and axial-vector, while the subscript refers to the spin J of 
the hadronic system. These spectral functions can be experimentally determined by 
measuring the invariant mass spectra of given hadronic modes and normalizing them to 
their respective branching ratios. The nonstrange spectral functions (fi/ai and ao) are 
defined and measured in [ (Lj, § , and the strange spectral functions read 



Ml 



B{t- -> V~/A~(S = -1, J = IK) 



6 \V US \ Sew 



1 diV, 



V/A 



and 



v${s)/a${s) 



V/A 



Ml 



ds 



B(t~ - 


-> e z/ e z/ T ) 




/ 2s 







(34) 



S(r- -> V-/A~(S — —1, J = 0)i/ T ) 



6 fKw 2 3£W 



Bfr- 



-i/ e z/ T ) 



1 rfAT, 



V/A 



1 



(35) 



where I V,, 



0.2218 ± 0.0016 

Sew — 



M T = (1776.91^) MeV/c 2 is taken from the 
BES measurement [ |36|, Sew — 1-0194 ± 0.0040 accounts for the electroweak radiative 
corrections [0, and dNy/A/ 'Ny/Ads denotes the normalized distribution of the invariant 
mass-squared of the corresponding vector /axial- vector decay channels V/A. The leptonic 
branching ratio is taken to be B(t~ — > e~v e v T ) = (17.794 ± 0.045)%, where the precision 
has been improved by applying lepton universality as described in Ref. [ 0]. The 

s) is expected to be provided by the single kaon pole with 

^P K 5{s-Ml). 

Because the spectral functions measure the transition probability to create hadrons 
with a mass yfs from the QCD vacuum, they are the natural input to QCD studies, 
allowing tests to be performed at a running mass scale less than or equal to M T . Their 
nonstrange vector part can be directly related to the low energy e + e~ annihilation cross 
section using CVC. In addition, the strange spectral function carries information on the 
strange quark mass. 



dominant contribution to 

a 0,K\ b , 
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Figure 9: Total hadronic vector and axial-vector spectral function (dots with error bars) 
from r decays into strange final states with its different contributions indicated. The errors 
include systematic uncertainties. The kaon pole is not shown in this plot. The shapes for 
the K~r] and Kmr(n > 3) modes are obtained from the Monte Carlo simulation. The 
parton model prediction is given by the solid straight line. 



The total strange spectral function is shown in Fig. |^. The spectra for the Kir and Kirn 
modes are obtained from the corrected mass spectra, taking into account the acceptance 
and bin migration corrections. For the other modes a phase space generator is used to 
simulate the corresponding mass distributions. The contributions from the respective 
channels are normalized to the corresponding branching ratios according to Eq. (51). 

The i^*(892)~ resonance stands out clearly on the low mass side, while within a large 
uncertainty (dominated by background subtraction in the (Kmr)^ modes) the higher 
energy part is consistent with the parton model expectation. 

The total strange spectral function is used in the following section, in order to construct 
spectral moments needed for the QCD analysis. 
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7 The strange quark mass 



A main goal of the present analysis is the determination of the strange quark mass 
m s . The method adopted follows the line of the recent ALEPH a s (M 2 ) determination 
from nonstrange hadronic r decays [ |2j and is based on a simultaneous fit of QCD 
parametrizations, including perturbative and nonperturbative components, to measured 
observables. A more detailed account of this method and of the phenomenological context 
can be found in Ref. 
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7.1 The strange hadronic decay ratio R T s 

|, the inclusive r decay ratio into strange hadronic 



As previously demonstrated in Ref. 
final states, 

Rr,t 



T(r — > hadrons 5= _ 1 v T ) 



rfr- 



e fpf-r 



(36) 



can be used due to its precise theoretical prediction [ [K| |TTJ in the framework of the 
Operator Product Expansion (OPE) [ to determine m s (s) at the scale s = M 2 . Since 
then it was shown [ that the perturbative expansion used for the massive term in Ref. [ 
TTJ was incorrect. After correction the series shows a problematic convergence behaviour [ 



12, 13, 14 



Following Ref. [ [T(| the theoretical prediction for the inclusive vector plus axial- vector 
hadronic decay rate is given by 



R T {M 2 T ) 



12vrS- 



EW 



Ml 

as 



Ml 




x 



1 + 2^ ImnW( S ) - 2^1 



s) 



(37) 



with the two-point correlation functions lT J ) = | V^^n^/y^ + \V^ 
hadronic final state of the spin J. The choice of the particular spin combination for the 
correlators taken in Eq. ([37]) is justified in Ref. [ [10| . Equation ([37]) can be decomposed 
as 



R 



T,(S) 



3\V\ 2 S EW 1 + 5 (0) + 5 (2 ~ mass } + 6(D) + s 'i 



EW 



(38) 



D=4,6, 



where V = V u d {V = V us ) for the nonstrange (strange) case is the corresponding CKM 
matrix element. The residual non-logarithmic electroweak correction 5' EW ~ 0.0010 [ [^J 
is neglected in the following. Throughout this analysis, all QCD observables are expressed 
in the MS renormalization scheme. 



7.1.1 Perturbative contributions 



The 5^ term in Eq. ( p8|) is the perturbative part of mass dimension D — 0, known to 
third order [ ^0f in the expansion with a s (M 2 ) = a s (M 2 )/ir, and studied in detail in Ref. [ 
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Next, the S^ 2 mass * ) term is the mass contribution of dimension D = 2, in practice 
only important for the strange quark mass and the subject of the present study. The 



fixed-order perturbation theory (FOPT) gives [ LI, 12 



8s 



(2— mass) 



' Ml 



16 

1 + — a s (M 2 ) + 46.00 a 2 s {M 2 T ) 
3 (i+o) 



+ ( 283.6 + ^>)ai(M!) 



(39) 



The third order coefficient x^ ^ occurs in the expansion of the massive J = 1 + 
correlator [p, [|2), while the J = correlator is known up to third order. Assuming the 



perturbative expansions to behave like a geometric series, one can estimate the unknown 
coefficient to be x^ +0 ^ ~ a4 1+ °' ) ( x 2 1+0 V a; i L+0 ' 1 ) = 165 ± 330, using as error twice the 
estimated contribution. Setting a s (M 2 ) = 0.334 [ @, Eq. © becomes 

5 ( 2 -mass) = - [1 + 0.57 + 0.52 + (0.49 ± 0.29)1 , (40) 

which converges badly. 

A so-called contour-improved FOPT [ [E| fyj analysis (FOPTci) for the dimension 
D = 2 contribution has been presented in Refs. [ [13], [14]]. It consists of a direct 
numerical evaluation of the contour integral derived from Eq. (|37D, using the solution of 
the renormalization group equation (RGE) to four loops [ £|5|, fJ6] as input for the running 
a s (s) and m s (s). This provides a resummation of all known higher order logarithmic 
integrals and was observed to improve the convergence of the massless perturbative series [ 
H, |43| . With the above value for a s , the contour-improved evaluation of the perturbative 



series reads 

m 2 ( M 2 ) 

^g-mass) = _ [0.97 + 0.49 + 0.38 + (0.34 ± 0.04)] , (41) 

with somewhat improved convergence compared to FOPT (flop. Therefore, the results of 
the following analysis are based on the FOPT CI approach. 

Independently of whether FOPT or FOPTci is used, the origin of the convergence 
problem is found in the J = component as defined in Eq. (|57D 

4 2 ~ mass) (</ = 0) = -8^|^ [0.41 + 0.32 + 0.32 + 0.37] , (42) 

while the J = 1 + part converges very well, 

4 2 " mass) (J= 1 + 0) = - 8 m 'fy [0.56 + 0.16 + 0.05 - (0.04 ± 0.04)] . (43) 

Following these observations, two methods are considered in the following in order to 
determine m s (M 2 ): 
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• in the inclusive method, the inclusive strange hadronic rate is considered and both 
J = 1 + and J = are included with their respective convergence behaviour taken 
into account in the theoretical uncertainties. 

• the '1+0' method singles out the well-behaved J = 1 + part by subtracting 
the experimentally determined J = longitudinal component from data. The 
measurement is then less inclusive and the sensitivity to m s is significantly reduced; 
however, the 62 perturbative expansion is under control and the corresponding 
theoretical uncertainty is reduced. 

Since the J = expansion is problematic, the results are given with the '1+0' method, 
while the inclusive approach is used as an insight into the handling of Eq. fl42|) . 



7.1.2 Nonperturbative contributions 



The last term in Eq. (38) represents the nonperturbative contribution which, using the 
OPE, can be written as a sum of powers of M®: 

dimO=D V 1V1 t) 

where the parameter fi separates the long-distance nonperturbative effects, absorbed 
into the vacuum expectation elements from the short-distance effects which are 

included in the Wilson coefficients C(s,/i) [ f47f . 

The dimension D = 4 operators have dynamical contributions from the gluon 
condensate {(a s /n)GG), the quark condensates m u (uu), m d (dd) and m s (ss), and the 
running quark masses to the fourth power. The contributions from dimension D = 6 
and D = 8 operators are rather complex and they are taken into account by introducing 
effective scale independent operators (Oq) and (Og) that are fitted to the data. 



7.2 Evidence for the effect of a massive strange quark 

From the result found from this analysis on the branching ratio for r decays into all 
strange hadronic final states (|2"7D, one obtains 

R T<S = 0.1610 ± 0.0066 , (45) 

using the value for B{r~ — > e~z/ e i/ r ) given in Section 6. 

The result (|45D can be readily compared with the QCD prediction for a massless 
strange quark neglecting the nonperturbative contributions, 

R { °1 = 0.1809 ±0.0036 , (46) 

where the quoted error mainly reflects the uncertainties on a s (M?) and V us . The two 
values differ by 2.7a, in the direction predicted for a non-zero m s value. 

This is evidence for the effect of a massive strange quark in a single observable which 
is essentially predicted by perturbative QCD, providing the basis for the determination 
of the running mass m s at the scale of the r mass. 
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7.3 Spectral moments 



As proposed in Ref. [[4l| and successfully applied in several a s (M^) analyses [ 0, [41], [51 
the spectral functions are used to construct the moments 



(47) 



with R T = Rt,(S) f° r the nonstrange and the strange cases, respectively. The theoretical 
prediction for the moments (|47|) follows the line described above, leading to expressions 
similar to Eq. (^) with nonperturbative contributions £ W ( D ) [|15|. Fitting the r decay rate 
and the spectral moments allows m s (M^) and the nonperturbative operators of dimension 
D = 6 and D = 8 to be simultaneously obtained. 

In order to reduce the theoretical uncertainties, the difference between nonstrange 
and strange spectral moments, properly normalized with their respective CKM matrix 
elements, is considered: 

A r = W^ R r,S=0 ~ T^T2 R T,S=-l , (48) 
| *ud | | * us | 

for which the massless perturbative contribution vanishes so that the theoretical prediction 
now reads (setting m u = m d = 0) 

A? = 3S EW (-5 k s l(2 - mass) + ~ 6kl(D) ) ■ ( 49 ) 

For the dimension D = 4 nonperturbative contribution, the flavour independent gluon 
condensate disappears in the difference, leaving only the quark condensate terms. The 
fitted 5 fc ^ 6 ' 8 ) contain the higher-order nonperturbative contributions to the difference of 
nonstrange and strange moments. Since the nonstrange nonperturbative contributions to 
the inclusive V + A hadronic decay ratio have been found to be very small [ 0] , one has 

^fc/(4,6,8) ^ _^M(4,6,8) 

In this analysis, the nonstrange spectral function obtained from Ref. [ [| has been 
appropriately scaled making use of the updated value of -R T ,s@ The experimental results 
of the spectral moments (fl"8|) and their correlations are given in Table || For the CKM 
matrix elements the values \V ud \ = 0.9751 ± 0.0004 and \V US \ = 0.2218 ± 0.0016 [ gU are 
used, while the errors are included in the theoretical uncertainties discussed in the next 
section. In the inclusive V + A nonstrange case the correlations between the hadronic 
decay ratio, obtained essentially using universality from the leptonic branching ratios and 
the r lifetime, and the spectral moments are negligible (see Ref. [0). On the contrary, 
all (k, I) moments used here suffer from large correlations due to the common input from 
the strange spectral function. This reduces the independent information used in the fit 
and thus generates strong correlations between the adjusted parameters. 

Figure [10] shows the weighted integrand of the lowest moment A°° (see Eq. (|4"7D , (pE8]) ) 
from the ALEPH data, as a function of the invariant mass-squared, and for which the 
expectation from massless perturbative QCD vanishes. 

4 i? Tj s=o is obtained from the difference R T — Rt,s, in which R T is directly derived from the leptonic 
branching ratio quoted in Section 6. 
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(k, I) 




(0,0) 


0.394 ±0.137 


(1,0) 


0.383 ±0.078 


(2,0) 


0.373 ± 0.054 


(1,1) 


0.010 ±0.029 


(1,2) 


0.006 ±0.015 



(k,l) 


(0,0) 


(1,0) 


(2,0) 


(1,1) 


(1,2) 


(0,0) 


1 


0.94 


0.83 


0.98 


0.91 


(1,0) 




1 


0.97 


0.87 


0.72 


(2,0) 






1 


0.72 


0.53 


(1,1) 








1 


0.95 


(1,2) 










1 



Table 8: Measured spectral moments A!f (left table) and their experimental correlations 
(right table). 
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Figure 10: Integrand of Eq. ffify for (k=0, 1=0), i.e., difference of the Cabibbo 
corrected nonstrange and strange invariant mass spectra. The contribution from massless 
perturbative QCD vanishes in the difference. 



7.4 Separation of J = 0, 1 components 

Although the J = component has not been thoroughly investigated in the data, it is 
clear that it is dominated by the single K pole. Other contributions can be identified: 
(i) off-shell vector K * resonances can generate a + component [ , however too small 
to have a significant effect in this analysis, and (ii) production of the scalar ^(1430) 
resonance may occur as discussed below. 

The _K"q(1430) state decays almost exclusively into with a branching ratio of 
(93 ± 10)% [ |18| and a fit to the K~tc° mass distribution as discussed in Section 4.1 
provides a branching ratio B(t — > /Tq(1430) v t ) = (0.0 ±2.5) x 10 -4 . It is therefore clear 
that + contributions are strongly suppressed (< 1.6 x 10~ 2 ) compared to the dominant 1~ 
production. This is expected from chiral symmetry breaking, with contributions reduced 
by ~ (m^/M r ) 4 as discussed in Ref. [ 0], not including the r kinematic factor which 
disfavours higher mass states. 
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No direct estimate of extra 0~ components beyond the single K can be made from 
the present data in the Ktttt modes. This is due to the lack of statistics and to the 
background subtraction, especially in the high mass region where some evidence for a 
broad pseudoscalar state at 1460 MeV/c 2 exists [ A contribution equal to the + one 
(actually an upper limit) is assumed in this case. 

The evaluation of the J = contributions to R Tt s obtained from integration of Eq. Q3"7|) 
using f K from Eq. @ gives 

R^ S (K) = - 0.00615 ± 0.00026 , (50) 
i4°s(other 0~,0 + ) = - 0.0015 ± 0.0015 , (51) 

leading to the result 

A°°(J = 0) = 0.155 ±0.031, (52) 

where the J = contributions in the nonstrange part have been computed to be —0.00076 
for the 7r pole (using the corresponding ALEPH branching ratio [|3~H1) and safely assumed 



to be negligible in the other modes. Similarly, the J = contributions are obtained for 
all (k, I) values and subtracted from the measured moments. 

7.5 Theoretical parameters and uncertainties 

The recent ALEPH measurement of the strong coupling constant [ |2[, a s (M^) = 
0.334 ± 0.007 cxp ± 0.021th , is used. The theoretical error is dominated by uncertainties 
due to the truncation of the perturbative series so that correlations between the above 
value and the moments fl4"8|) are negligible, with the latter dominated by the uncertainties 
on the strange part. 

The dependence on the renormalization scheme (MS) used is already included in the 
theoretical error on a s (M^) and thus not further added in order to avoid double counting. 
The renormalization scale (/i) dependence of the D = 2 prediction has been studied 



in Ref. [ [L3[. When varying fi away from M T , additional logarithms enter the series and 
due to its truncation at finite order, a residual dependence on \i is left. To estimate the 
associated uncertainty, fi is varied from 1.3 to 3.6 GeV/c 2 . 

The uncertainty originating from the truncation of the perturbative series for the 
dimension D = 2 term is handled as follows: the estimated error is taken as the full size 
of the last term retained in the expansion, namely of order for the J = 1 + series. 
Some redundancy is expected between the uncertainties estimated for the renormalization 
scale dependence and the truncation of the perturbative series. However, both errors are 
conservatively kept and added in quadrature. 

For the dimension D = 4 nonperturbative contribution the products of quark 
condensates and quark masses are taken from the PCAC relations with some correction 
for the strange quark [ No errors are introduced for the higher dimensional operators 
since they are fixed experimentally. Finally, a small uncertainty is included from the error 
on Sew- 

Table [5] gives the above theoretical uncertainties on m s (M^). The main error 
contribution stems from the |V^ S | uncertainty 
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truncation 




11 


R-scale \x 


(1.3 - 3.6) GeV/c 2 


6 


m s (ss) 


-(1.63 ±0.29) x 1(T 3 GeV 4 


1 


Sew 


1.0194 ±0.0040 


1 


\v us \ 


0.2218 ±0.0016 
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Total errors 




26 



Table 9: Theoretical uncertainties on m s (M 2 ) . The uncertainty from the truncation of the 
perturbative series for the D = 2 mass term is estimated from adding or subtracting the 
value of the last retained term in the QCD expansion. The unequal positive and negative 
errors on m s are averaged. 





m s (M 2 ) 6® 


s® 


m s (M 2 ) 


1 0.74 


-0.92 




1 


-0.83 


#8) 




1 



Table 10: Correlation matrix according to the '1+0' fit results. 



7.6 Results 

The fit minimizes the y 2 of the differences between measured and fitted quantities 
contracted with the inverse of the sum of the experimental and theoretical covariance 
matrices. Due to the large correlations and the statistical limitation of the spectral 
functions at the end of the r phase space, higher moments / > 2 do not add significant 
information to the fit. The results for m s (M 2 ) and the nonperturbative contributions to 
A oo(i+o) for the t 1+0 > method are 

m s (M 2 T ) = (176^1^ ± 8 fit ± Hj=o) MeV/c 2 , (53) 
6® = 0.039 ± 0.016 cxp ± 0.014th ± 0.004 fit , (54) 
5® = - 0.021 ± 0.014 cxp ± 0.008 t h ± 0.003 fit , (55) 

with a x 2 /ndf of 0.2/2. 

The errors are separated according to their experimental, theoretical, fit and spin 
separation origins. The fit error stems from an intrinsic bias in the x 2 minimization, due 
to the large correlations of the input observables, and includes the total difference between 
the fit results with and without correlations, while the results without correlations are 



given as central values, as discussed in Ref . [ |53[ . The correlation matrix corresponding 



to the fit results is given in Table |10. 
The result fl53|) 

m s (M 2 ) = (176 tf 7 ) MeV/c 2 (56) 
can be evolved to the scale of 1 GeV using the four- loop RGE 7-function [ 46]], yielding 

m s (l GeV 2 ) = (234+^) MeV/c 2 . (57) 
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The fitted m s value corresponds to a contribution 5 S = —(0.058 ± 0.027) to R T s . 
No operator of dimension D = 4 has been fitted except the small quartic strange mass 
corrections, for a total contribution 8*<p = —(0.003 ± 0.001), dominated by the strange 
quark condensate; the error given is almost entirely of theoretical origin. The D = 6 
and D = 8 strange contributions are 4 6) = -0.038 ± 0.022 and 4 §) = 0-°20 ± 0.016. 
They are fairly large compared to the nonstrange case, <5 (6) = 0.001 ± 0.004 and 
= -0.001 ±0.001 [g. 

The above results are obtained from the J = 1 + piece which was shown in Section 
7.1.1 to have a satisfactory convergence in QCD. The inclusive method has a better 
experimental definition, but suffers from the bad convergence properties of the J = part 
for which a prescription must be specified. A reasonable rule to handle an asymptotic 
series is to truncate it where the terms reach a minimum and then assign as an uncertainty 
the full amount of the last term retained. With this prescription, the J = expansion is 
stopped after the 0(a s ) term and, going through the same steps as before, one gets 

m s (M 2 T ) = (149^1^ ±6 fit ) MeV/c 2 , (58) 

where this time the theoretical uncertainty is dominated by the truncation of the 
perturbative series. The fitted nonperturbative parts are = 0.031 ± 0.015 and 
<J(8) = -0.018 ± 0.011. The results (0) and fl58|) are consistent within the uncorrelated 
part of their errors with a x 2 / n df of 1.4/1. 

It may be interesting to come back to the discarded J = part in order to get some 
information regarding the handling of its ill-behaved perturbative expansion, Eq. ([1*3]). 
Truncating the series at the minimum provides a result consistent with (|53|), while 
keeping the two more known terms destroys the consistency. This observation supports 
the prescription used in the inclusive method. However, as stated above, the safer 
determination with the '1+0' method is preferred and given as the final result of this 
analysis. 

As shown for the nonstrange case in Ref. [ 0, one can simulate the physics of a 
hypothetical r lepton with mass s/so smaller than M T by replacing M 2 everywhere 
in Eq. fl37|) by s . Under the assumption of quark-hadron duality, the evaluation of 
the observables as a function of so constitutes a test of the OPE approach adopted 
here, since the energy dependence of the theoretical predictions is determined once the 
parameters are fixed. Figure ITT] shows the running observable A°°( 1+ °) (sq) compared to the 
corresponding theoretical predictions and the fitted parameters in Eq. (|53],[53|[55[) . Despite 
the expected breakdown of the perturbative approach at lower scales, it is noteworthy 
that the prediction agrees with data at such a low mass scale within the theoretical 
uncertainties used in the analysis. 



7.7 Comparison with other determinations of m s 

Other determinations of m s have been obtained by analyses of the divergence of the vector 
and axial- vector current two-point function correlators [ [54] 
|B"3"| , |B4[ ]. The phenomenological information on the associated scalar and pseudoscalar 
spectral functions is reconstructed from phase-shift resonance analyses which are yet 
incomplete over the considered mass range and need to be supplemented by other assumed 



56, 57, m pa eq, bt 
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Figure 11: The observable A°°( 1+ °) (so) as a function of the "t -mass" -squared so. The 
curve is plotted as a one-standard deviation error band to emphasize its strong point-to- 
point correlations. Also shown is the theoretical prediction for the fit parameters given in 
Eq. gS3J5|). 



ingredients, in particular the description of the continuum, thus introducing systematic 
effects. 

Another approach [ considers the difference between isovector and hypercharge 
vector current correlators as related to the 1 = 1 and 1 = spectral functions accessible 
in e + e~ annihilation into hadrons at low energy. A recent reanalysis [ points out the 
possibility of large corrections from isospin breaking, leading to significant deviations for 
the extracted m s value. 

Finally, lattice QCD calculations of m s are available [ |67|, |68], |69], [70||. Their values 
show quite a large spread. 

The present determination of m s is directly compared with the other derivations in 



Fig. [12] where all values are given at a mass scale of 1 GeV. Qualitative agreement is 
observed, the value from the present analysis being rather on the high side of the range of 
previous determinations. The precision of the result is unfortunately still limited. As an 
example, the Standard Model prediction for the CP- violation parameter e'/e is strongly 
dependent on the value of m s and the present determination translates into a range of 
e'/e values spanning a factor of five [ 0], but favouring "low" values for this ratio. 



8 Conclusion 

All ALEPH measurements on r decays with kaons are summarized to provide an overall 
review and a comprehensive study of these decays. 

Comparing the decays r _ — ► K~v T and K~ —>■ n~v T shows agreement with fi — r 
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Figure 12: The ALEPH determination of m s (l GeV 2 ) compared to the results of other 
approaches. Details are given in the text (SR = sum rules). The references are listed in 
the order of the results from top to bottom. 

universality within 2% uncertainty. 

The investigation of mass spectra shows that the (Kn)~ decay mode is dominated by 
if* (892) ~. A if* (1410) contribution is extracted from a fit to the Kn mass spectrum 
through its interference with the dominant if* (892). 

The (Kmr)^ decay modes are observed to proceed through the ifi(1270), ifi(1400) 
and if* (1410) resonances. The relative contributions are extracted by fitting the 7T7T 
and Kixix mass spectra, separating the vector and axial-vector contributions. The vector 
current contribution to the Kirn modes is found to be (23t 2 j)%. 

According to CVC, the vector spectral function for the {KKix)~ modes is connected 
to low energy e + e~ annihilation, allowing one to conclude that (94lg)% of this decay 
mode proceeds through the axial- vector current. Because the latter is dominated by the 
a\ resonance, a branching ratio B(ai — > K*K) = (2.6 ± 0.3)% is obtained. 

Tests of isospin symmetry for the (if7r)~, (Knir)' and (KK-ir)~ modes are performed. 
A good consistency with the expectation from isospin invariance is observed, providing 
useful information on the decay dynamics in these final states. 

The branching ratio for r decay into all strange final states is determined to be 
B(t~ — >• X~(S = -l)z/r) = (28.7 ± 1.2) x 10~ 3 and the total strange spectral function is 
derived from the corresponding mass spectra. 

Using the total rates, as well as moments of the respective spectral functions, through 
a combination of strange and nonstrange parts which cancel the dominant massless 
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perturbative QCD contribution, a global fit is performed providing a determination of 
the strange nonperturbative contributions and of the s quark mass. The value obtained 
is m s (M 2 ) = (176 If?) MeV/c 2 , which is evolved to 1 GeV to yield m s (l GeV 2 ) = 
(234 +7g) MeV/c 2 , in agreement with other determinations. 
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